Anthropogenic contamination sources drive differences in antimicrobial-resistant Escherichia coli in three urban lakes.

Antimicrobial resistance (AMR) is an ever-present threat to the treatment of infectious diseases. However, the potential relevance of this phenomenon in environmental reservoirs still raises many questions. Detection of antimicrobial-resistant bacteria in the environment is a critical aspect for understanding the prevalence of resistance outside of clinical settings, as detection in the environment indicates that resistance is likely already widespread. We isolated antimicrobial-resistant Escherichia coli from three urban waterbodies over a 15-month time series, determined their antimicrobial susceptibilities, investigated their population structure, and identified genetic determinants of resistance. We found that E. coli populations at each site were composed of different dominant phylotypes and showed distinct patterns of antimicrobial and multidrug resistance, despite close geographic proximity. Many strains that were genome-sequenced belonged to sequence types of international concern, particularly the ST131 clonal complex. We found widespread resistance to clinically important antimicrobials such as amoxicillin, cefotaxime, and ciprofloxacin, but found that all strains were susceptible to amikacin and the last-line antimicrobials meropenem and fosfomycin. Resistance was most often due to acquirable antimicrobial resistance genes, while chromosomal mutations in gyrA, parC, and parE conferred resistance to quinolones. Whole-genome analysis of a subset of strains further revealed the diversity of the population of E. coli present, with a wide array of AMR and virulence genes identified, many of which were present on the chromosome, including blaCTX-M. Finally, we determined that environmental persistence, transmission between sites, most likely mediated by wild birds, and transfer of mobile genetic elements likely contributed significantly to the patterns observed.IMPORTANCEA One Health perspective is crucial to understand the extent of antimicrobial resistance (AMR) globally, and investigation of AMR in the environment has been increasing in recent years. However, most studies have focused on waterways that are directly polluted by sewage, industrial manufacturing, or agricultural activities. Therefore, there remains a lack of knowledge about more natural, less overtly impacted environments. Through phenotypic and genotypic investigation of AMR in Escherichia coli, this study adds to our understanding of the extent and patterns of resistance in these types of environments, including over a time series, and showed that complex biotic and abiotic factors contribute to the patterns observed. Our study further emphasizes the importance of incorporating the surveillance of microbes in freshwater environments in order to better comprehend potential risks for both human and animal health and how the environment may serve as a sentinel for potential future clinical infections.

Recurring Trans-Atlantic Incursion of Clade 2.3.4.4b H5N1 Viruses by Long Distance Migratory Birds from Northern Europe to Canada in 2022/2023.

In December 2022 and January 2023, we isolated clade 2.3.4.4b H5N1 high-pathogenicity avian influenza (HPAI) viruses from six American crows (Corvus brachyrhynchos) from Prince Edward Island and a red fox (Vulpes vulpes) from Newfoundland, Canada. Using full-genome sequencing and phylogenetic analysis, these viruses were found to fall into two distinct phylogenetic clusters: one group containing H5N1 viruses that had been circulating in North and South America since late 2021, and the other one containing European H5N1 viruses reported in late 2022. The transatlantic re-introduction for the second time by pelagic/Icelandic bird migration via the same route used during the 2021 incursion of Eurasian origin H5N1 viruses into North America demonstrates that migratory birds continue to be the driving force for transcontinental dissemination of the virus. This new detection further demonstrates the continual long-term threat of H5N1 viruses for poultry and mammals and the subsequent impact on various wild bird populations wherever these viruses emerge. The continual emergence of clade 2.3.4.4b H5Nx viruses requires vigilant surveillance in wild birds, particularly in areas of the Americas, which lie within the migratory corridors for long-distance migratory birds originating from Europe and Asia. Although H5Nx viruses have been detected at higher rates in North America since 2021, a bidirectional flow of H5Nx genes of American origin viruses to Europe has never been reported. In the future, coordinated and systematic surveillance programs for HPAI viruses need to be launched between European and North American agencies.

Extracellular Polysaccharide Receptor and Receptor-Binding Proteins of the Rhodobacter capsulatus Bacteriophage-like Gene Transfer Agent RcGTA.

A variety of prokaryotes produce a bacteriophage-like gene transfer agent (GTA), and the alphaproteobacterial Rhodobacter capsulatus RcGTA is a model GTA. Some environmental isolates of R. capsulatus lack the ability to acquire genes transferred by the RcGTA (recipient capability). In this work, we investigated the reason why R. capsulatus strain 37b4 lacks recipient capability. The RcGTA head spike fiber and tail fiber proteins have been proposed to bind extracellular oligosaccharide receptors, and strain 37b4 lacks a capsular polysaccharide (CPS). The reason why strain 37b4 lacks a CPS was unknown, as was whether the provision of a CPS to 37b4 would result in recipient capability. To address these questions, we sequenced and annotated the strain 37b4 genome and used BLAST interrogations of this genome sequence to search for homologs of genes known to be needed for R. capsulatus recipient capability. We also created a cosmid-borne genome library from a wild-type strain, mobilized the library into 37b4, and used the cosmid-complemented strain 37b4 to identify genes needed for a gain of function, allowing for the acquisition of RcGTA-borne genes. The relative presence of CPS around a wild-type strain, 37b4, and cosmid-complemented 37b4 cells was visualized using light microscopy of stained cells. Fluorescently tagged head spike fiber and tail fiber proteins of the RcGTA particle were created and used to measure the relative binding to wild-type and 37b4 cells. We found that strain 37b4 lacks recipient capability because of an inability to bind RcGTA; the reason it is incapable of binding is that it lacks CPS, and the absence of CPS is due to the absence of genes previously shown to be needed for CPS production in another strain. In addition to the head spike fiber, we found that the tail fiber protein also binds to the CPS.

Distinct epidemiological profiles of porcine circovirus 3 and fox circovirus in Canadian foxes (Vulpes spp.).

Circoviruses (genus Circovirus, family Circoviridae) are ssDNA viruses that infect mammals, and they sometimes can transmit among different species. We investigated the distribution and diversity of porcine circovirus 3 (PCV-3, species Porcine circovirus 3) and fox circovirus (species Canine circovirus 1) in different populations of foxes (Vulpes spp.) inhabiting the Canadian province of Newfoundland and Labrador to compare their epidemiological profiles. Of the 210 samples tested in this study 9 were positive for PCV-3 and 99 were positive for fox circovirus. Eight foxes were PCV-3-positive (8/128, 6.3%) and the virus was only found in the most human-populated areas. The PCV-3 positivity rate was significantly higher in stool (7/180, 8.8%) than in spleen (2/120, 1.7%: p < 0.05). Phylogenetic analyses showed that sequences from different animals were unrelated to each other. Fox circovirus was identified in 66 animals (51.6%) and positivity rates were the highest in the least human-populated areas. There were no significant differences between positivity rates in stool (32/80, 40.0%), spleen (59/120, 49.2%), or lymph nodes (8/10, 80.0%). Among 54 positive animals for which both spleen and stool samples were available, 25 (46.3%) had detectable virus in both samples. All fox circovirus sequences recovered in this study formed a monophyletic clade, and no geographic segregation of study strains was observed. The high prevalence and high genetic diversity observed for fox circovirus implies that the virus has been circulating in this population for a long time. PCV-3 cases were consistent with sporadic infections from multiple sources, possibly related to scavenging behavior and consumption of meat by-products and human waste, while fox circovirus was endemic, indicating that foxes are likely the maintenance host for this virus. To the best of our knowledge this is the first study demonstrating the presence of fox circovirus in North America and to show that PCV-3 can be detected in foxes. Future studies should evaluate the pathogenic potential of these viruses for wildlife.

A novel calicivirus discovered in trumpeter swans (Cygnus buccinator) expands the richness of known avian caliciviruses.

Caliciviruses are ssRNA viruses that can infect a wide range of hosts, including birds. While several avian caliciviruses have been discovered, their taxonomy and host distribution are largely unknown. We molecularly characterized a novel calicivirus (trumpeter swan calicivirus: TruSCV) in trumpeter swans over-wintering in south-west British Columbia, Canada. The positivity rate was 20.3% (14/69) and there were no significant differences in infection rates between males (5/34, 14.7%) and females (9/35, 25.7%) or among considered age groups (juveniles: 4/14, 28.6%; sub-adults: 1/9, 11.1%; adults: 9/46, 19.6%). Twelve infected swans died of lead poisoning, one because of starvation, and one from physical injuries. TruSCV complete genome possessed the typical organization and protein motifs of caliciviruses and a type 2 IRES and its closest relative was a virus circulating in Australian ducks. Phylogenetic analyses showed the existence of 34 different but monophyletic avian caliciviruses. These viruses, while having conserved genomic organization and protein motifs, possess different IRES types and group in several divergent clades, with only two of them corresponding to currently defined genera, highlighting the need for epidemiological investigations and systematic analyses to better define their taxonomy. Follow-up studies are needed to elucidate the diversity, distribution, and pathogenic potential of TruSCV.

Patterns of abundance, chromosomal localization, and domain organization among c-di-GMP-metabolizing genes revealed by comparative genomics of five alphaproteobacterial orders.

BACKGROUND: Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is a bacterial second messenger that affects diverse processes in different bacteria, including the cell cycle, motility, and biofilm formation. Its cellular levels are controlled by the opposing activities of two types of enzymes, with synthesis by diguanylate cyclases containing a GGDEF domain and degradation by phosphodiesterases containing either an HD-GYP or an EAL domain. These enzymes are ubiquitous in bacteria with up to 50 encoded in some genomes, the specific functions of which are mostly unknown. RESULTS: We used comparative analyses to identify genomic patterns among genes encoding proteins with GGDEF, EAL, and HD-GYP domains in five orders of the class Alphaproteobacteria. GGDEF-containing sequences and GGDEF-EAL hybrids were the most abundant and had the highest diversity of co-occurring auxiliary domains while EAL and HD-GYP containing sequences were less abundant and less diverse with respect to auxiliary domains. There were striking patterns in the chromosomal localizations of the genes found in two of the orders. The Rhodobacterales' EAL-encoding genes and Rhizobiales' GGDEF-EAL-encoding genes showed opposing patterns of distribution compared to the GGDEF-encoding genes. In the Rhodobacterales, the GGDEF-encoding genes showed a tri-modal distribution with peaks mid-way between the origin (ori) and terminus (ter) of replication and at ter while the EAL-encoding genes peaked near ori. The patterns were more complex in the Rhizobiales, but the GGDEF-encoding genes were biased for localization near ter. CONCLUSIONS: The observed patterns in the chromosomal localizations of these genes suggest a coupling of synthesis and hydrolysis of c-di-GMP with the cell cycle. Moreover, the higher proportions and diversities of auxiliary domains associated with GGDEF domains and GGDEF-EAL hybrids compared to EAL or HD-GYP domains could indicate that more stimuli affect synthesis compared to hydrolysis of c-di-GMP.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria.

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

Virus-derived gene transfer agents benefit host cells by providing templates for DNA repair.

The potential benefits of producing gene transfer agents (GTAs) have long been speculated. A new study in PLOS Biology shows that DNA transfer by these phage-like elements allows cells with DNA damage to perform DNA repair and survive.

Formal recognition and classification of gene transfer agents as viriforms.

Morphological and genetic features strongly suggest that gene transfer agents (GTAs) are caudoviricete-derived entities that have evolved in concert with cellular genomes to such a degree that they should not be considered viruses. Indeed, GTA particles resemble caudoviricete virions, but, in contrast to caudoviricetes (or any viruses), GTAs can encapsidate at best only part of their own genomes, are induced solely in small subpopulations of prokaryotic host cells, and are transmitted vertically as part of cellular genomes during replication and division. Therefore, the lifecycles of GTAs are analogous to virus-derived entities found in the parasitoid wasps, which have recently been recognized as non-virus entities and therefore reclassified as viriforms. We evaluated three distinct, independently exapted GTA groups, for which the genetic basis for GTA particle production has been established. Based on the evidence, we outline a classification scheme for these viriforms.

A new perspective on the evolution and diversity of the genus Amdoparvovirus (family Parvoviridae) through genetic characterization, structural homology modeling, and phylogenetics.

Amdoparvoviruses (genus Amdoparvovirus, family Parvoviridae) are primarily viruses of carnivorans, but recent studies have indicated that their host range might also extend to rodents and chiropterans. While their classification is based on the full sequence of the major nonstructural protein (NS1), several studies investigating amdoparvoviral diversity have been focused on partial sequences, leading to difficulties in accurately determining species demarcations and leaving several viruses unclassified. In this study, while reporting the complete genomic sequence of a novel amdoparvovirus identified in an American mink (British Columbia amdoparvovirus, BCAV), we studied the phylogenetic relationships of all amdoparvovirus-related sequences and provide a comprehensive reevaluation of their diversity and evolution. After excluding recombinant sequences, phylogenetic and pairwise sequence identity analyses allowed us to define fourteen different viruses, including the five currently classified species, BCAV, and four additional viruses that fulfill the International Committee on Taxonomy of Viruses criteria to be classified as species. We show that the group of viruses historically known as Aleutian mink disease virus (species Carnivore amdoparvovirus 1) should be considered as a cluster of at least four separate viral species that have been co-circulating in mink farms, facilitating the occurrence of inter-species recombination. Genome organization, splicing donor and acceptor sites, and protein sequence motifs were surprisingly conserved within the genus. The sequence of the major capsid protein virus protein 2 (VP2) was significantly more conserved between and within species compared to NS1, a phenomenon possibly linked to antibody-dependent enhancement (ADE). Homology models suggest a remarkably high degree of conservation of the spikes located near the icosahedral threefold axis of the capsid, comprising the surface region associated with ADE. A surprisingly high number of divergent amino acid positions were found in the luminal threefold and twofold axes of the capsid, regions of hitherto unknown function. We emphasize the importance of complete genome analyses and, given the marked phylogenetic inconsistencies across the genome, advise to obtain the complete coding sequences of divergent strains. Further studies on amdoparvovirus biology and structure as well as epidemiological and virus discovery investigations are required to better characterize the ecology and evolution of this important group of viruses.

Distribution and diversity of dog parvoviruses in wild, free-roaming and domestic canids of Newfoundland and Labrador, Canada.

Some parvoviruses of carnivorans can infect multiple host species. Since many canine parvoviruses were only discovered recently, their host-range is still unexplored. We examined the host distribution and diversity of five dog parvoviruses in four canine populations from Newfoundland and Labrador, Canada, and investigated the potential for these viruses to cross the species barriers. Canine bocavirus 2 (CBoV-2) and the minute virus of canines were detected in stool from free-roaming dogs from Labrador (5/48 [10.4%] and 3/48 [6.3%], respectively) and two different CBoV-2 variants were identified. Canine bufavirus was identified in stool from free-roaming dogs (1/48, 2.1%) and foxes (3/80, 3.8%) from Labrador, but two different variants were observed in the two host species. The variant found in foxes was highly divergent from previously identified strains. Two cachavirus 1 variants, genetically similar to those circulating in other Canadian wildlife, were found in spleens from Newfoundland coyotes (3/87, 3.5%). Canine parvovirus type 2 (CPV-2) was found in stool from free-roaming dogs from Labrador (2/48, 4.2%) and in spleens from Newfoundland coyotes (3/87, 3.5%). Comparing CPV-2 sequences from these hosts to those retrieved from local symptomatic domestic dogs revealed the presence of a highly heterogeneous viral population as detected strains belonged to five different clades. The close relationship between CPV-2a strains from a dog and a coyote suggests the occurrence of viral transfer between wild and domestic canids. The identification of highly related strains with a similar molecular signature characteristic of older CPV-2 strains in free-roaming and domestic dogs suggests a probable common ancestry and that older CPV-2 strains, which have not been identified in dogs since the 1990s, persist in this part of Canada. Follow-up studies should evaluate samples from a larger number of animals and host species to extensively investigate the possible occurrence of cross-species transmission for recently discovered parvoviruses.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA.

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Natural disease and evolution of an Amdoparvovirus endemic in striped skunks (Mephitis mephitis).

Striped skunks (Mephitis mephitis) densely populate the human-animal interface of suburbia throughout North America. Skunks share that habitat with numerous related mesocarnivores, where increased contact, competition for shared food and water sources and other stressors contribute to increased exposure and susceptibility to viral infection. The recently identified skunk amdoparvovirus (SKAV) has been detected at high prevalence in skunks and occasionally in mink, but its distribution in North America is unknown. To understand the impact of SKAV in striped skunks and the risk posed to related species, we investigated the geographic distribution of SKAV, analysed its genetic diversity and evolutionary dynamics and evaluated viral distribution in tissues of infected animals to identify possible mechanisms of transmission. SKAV was detected in 72.5% (37/51) skunks and was present at high rates at all locations tested across North America. Analysis of the complete genomic sequence of 29 strains showed a clear geographic segregation, frequent recombination and marked differences in the evolutionary dynamics of the major structural (VP2) and non-structural (NS1) proteins. NS1 was characterized by a higher variability and a higher percentage of positively selected codons. This could indicate that antibody-mediated enhancement of infection occurs in SKAV, an infection strategy that may be conserved across amdoparvoviruses. Finally, in situ hybridization revealed virus in epithelium of the gastrointestinal tract, urinary tract and skin, indicating that viral transmission could occur via oronasal, faecal and/or urinary secretions, as well as from skin and hair. The endemicity of SKAV over large geographic distances and its high genetic diversity suggest a long-term virus-host association. Persistent shedding and high environmental stability likely contribute to efficient viral spread, simultaneously offering opportunities for cross-species transmission with consequent risk to sympatric species, including domestic animals and wildlife.

Co-circulation of five species of dog parvoviruses and canine adenovirus type 1 among gray wolves (Canis lupus) in northern Canada.

Several viruses can infect wild carnivores but their impact on wildlife health is poorly understood. We investigated the presence, diversity and distribution of various DNA viruses in 303 wolves inhabiting a vast area of the Northwest Territories, Canada, over a period of 13 years. We found evidence for the presence of canine bufavirus (CBuV, 42.6%), canine parvovirus 2 (CPV-2, 34.0%), canine bocavirus 2 (CBoV-2, 5.0%), cachavirus (CachaV-1, 2.6%), canine adenovirus 1 (CAdV-1, 1%) and minute virus of canines (MVC, 0.3%). To our knowledge, this is the first detection of CBoV-2, MVC and CachV-1 in wild animals. We also demonstrate that CBuV and CachaV-1 were already circulating among wild animals at least 11 and 10 years, respectively, before their discoveries. Although CBuV prevalence was higher, CPV-2 was the most prevalent virus among juveniles, while CBuV infection was associated with poor nutrition conditions. Even if its prevalence was low, CachaV-1 had the highest multiple infection rate (87.5%). CadV-1 and MVC sequences were highly identical to reference strains, but we observed a high diversity among the other viruses and detected three new variants. One CPV-2 variant and one CBuV variant were endemic since the beginning of the 2000s in the entire investigated region, whereas one CBuV variant and two CBoV-2 variants were found in a more restricted area over multiple years and CachaV-1 was found only in one region. Two CPV-2 variants and one CachaV-1 variant were observed only once, indicating sporadic introductions or limited circulation. Different patterns of endemicity might indicate that viruses were introduced in the wolf population at different timepoints and that mixing between wolf packs may not be constant. Different epidemiological behaviors depend on viral factors like infectivity, transmission routes, pathogenicity and tissue-tropism, and on host factors like proximity to densely populated areas, carnivory and pack density and mixing.

Transmission from vaccinated individuals in a large SARS-CoV-2 Delta variant outbreak.

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response.

Purification of Functional Gene Transfer Agents Using Two-Step Preparative Monolithic Chromatography.

Background: Gene transfer agents (GTAs) are phage-like particles that transfer cellular genomic DNA between cells. A hurdle faced in studying GTA function and interactions with cells is the difficulty in obtaining pure and functional GTAs from cultures. Materials and Methods: We used a novel two-step method for purification of GTAs from R. capsulatus by monolithic chromatography. Results: Our efficient and simple process had advantages compared to previous approaches. The purified GTAs retained gene transfer activity and the packaged DNA could be used for further studies. Conclusions: This method is applicable to GTAs produced by other species and small phages, and could be useful for therapeutic applications.

Novel Divergent Members of the Kitrinoviricota Discovered through Metagenomics in the Intestinal Contents of Red-Backed Voles (Clethrionomys gapperi).

Metagenomic methods are powerful tools to investigate viral diversity in biological or environmental samples and to identify previously unknown viruses. We used RNA metagenomics to identify, in the gut of red-backed voles, the nearly complete genomes of two novel members of the Kitrinoviricota, a phylum including viruses with positive-sense ssRNA genomes encoding an RNA-directed RNA polymerase. The genome of a novel member of the Tombusviridae presented four open reading frames (ORFs); a -1 frameshift is potentially involved in generating the viral replicase. This sequence was part of a phylogenetic clade that did not include any officially classified species. The second genome presented a large ORF coding for a viral polyprotein containing the typical protein domains common to flexiviruses. The sequence clustered with currently known members of the Deltaflexiviridae. Both viruses appear to represent the first members of novel species in yet undefined genera. The identified viruses likely originated from the vole diet as members of the two viral families are known to infect plants and fungi, respectively. Investigating public databases demonstrated that a much higher richness than currently recognized exists for these two viral families, highlighting the need to update taxonomy systems and possibly also include genomes identified through metagenomics.

Prolactin and prolactin receptor expression in the HPG axis and crop during parental care in both sexes of a biparental bird (Columba livia).

During breeding, multiple circulating hormones, including prolactin, facilitate reproductive transitions in species that exhibit parental care. Prolactin underlies parental behaviors and related physiological changes across many vertebrates, including birds and mammals. While circulating prolactin levels often fluctuate across breeding, less is known about how relevant target tissues vary in their prolactin responsiveness via prolactin receptor (PRLR) expression. Recent studies have also investigated prolactin (PRL) gene expression outside of the pituitary (i.e., extra-pituitary PRL), but how PRL gene expression varies during parental care in non-pituitary tissue (e.g., hypothalamus, gonads) remains largely unknown. Further, it is unclear if and how tissue-specific PRL and PRLR vary between the sexes during biparental care. To address this, we measured PRL and PRLR gene expression in tissues relevant to parental care, the endocrine reproductive hypothalamic-pituitary- gonadal (HPG) axis and the crop (a tissue with a similar function as the mammalian mammary gland), across various reproductive stages in both sexes of a biparental bird, the rock dove (Columba livia). We also assessed how these genes responded to changes in offspring presence by adding chicks mid-incubation, simulating an early hatch when prolactin levels were still moderately low. We found that pituitary PRL expression showed similar increases as plasma prolactin levels, and detected extra-pituitary PRL in the hypothalamus, gonads and crop. Hypothalamic and gonadal PRLR expression also changed as birds began incubation. Crop PRLR expression correlated with plasma prolactin, peaking when chicks hatched. In response to replacing eggs with a novel chick mid-incubation, hypothalamic and gonadal PRL and PRLR gene expression differed significantly compared to mid-incubation controls, even when plasma prolactin levels did not differ. We also found sex differences in PRL and PRLR that suggest gene expression may allow males to compensate for lower levels in prolactin by upregulating PRLR in all tissues. Overall, this study advances our understanding of how tissue-specific changes in responsiveness to parental hormones may differ across key reproductive transitions, in response to offspring cues, and between the sexes.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant quickly rose to dominance in mid-2021, displacing other variants, including Alpha. Studies using data from the United Kingdom and India estimated that Delta was 40-80% more transmissible than Alpha, allowing Delta to become the globally dominant variant. However, it was unclear if the ostensible difference in relative transmissibility was due mostly to innate properties of Delta's infectiousness or differences in the study populations. To investigate, we formed a partnership with SARS-CoV-2 genomic surveillance programs from all six New England US states. By comparing logistic growth rates, we found that Delta emerged 37-163% faster than Alpha in early 2021 (37% Massachusetts, 75% New Hampshire, 95% Maine, 98% Rhode Island, 151% Connecticut, and 163% Vermont). We next computed variant-specific effective reproductive numbers and estimated that Delta was 58-120% more transmissible than Alpha across New England (58% New Hampshire, 68% Massachusetts, 76% Connecticut, 85% Rhode Island, 98% Maine, and 120% Vermont). Finally, using RT-PCR data, we estimated that Delta infections generate on average ∼6 times more viral RNA copies per mL than Alpha infections. Overall, our evidence indicates that Delta's enhanced transmissibility could be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on the underlying immunity and behavior of distinct populations.

Evidence of transmission from fully vaccinated individuals in a large outbreak of the SARS-CoV-2 Delta variant in Provincetown, Massachusetts.

Multiple summer events, including large indoor gatherings, in Provincetown, Massachusetts (MA), in July 2021 contributed to an outbreak of over one thousand COVID-19 cases among residents and visitors. Most cases were fully vaccinated, many of whom were also symptomatic, prompting a comprehensive public health response, motivating changes to national masking recommendations, and raising questions about infection and transmission among vaccinated individuals. To characterize the outbreak and the viral population underlying it, we combined genomic and epidemiological data from 467 individuals, including 40% of known outbreak-associated cases. The Delta variant accounted for 99% of sequenced outbreak-associated cases. Phylogenetic analysis suggests over 40 sources of Delta in the dataset, with one responsible for a single cluster containing 83% of outbreak-associated genomes. This cluster was likely not the result of extensive spread at a single site, but rather transmission from a common source across multiple settings over a short time. Genomic and epidemiological data combined provide strong support for 25 transmission events from, including many between, fully vaccinated individuals; genomic data alone provides evidence for an additional 64. Together, genomic epidemiology provides a high-resolution picture of the Provincetown outbreak, revealing multiple cases of transmission of Delta from fully vaccinated individuals. However, despite its magnitude, the outbreak was restricted in its onward impact in MA and the US, likely due to high vaccination rates and a robust public health response.